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Abstract 

We analyse the "difFractive" events reported by the HERA groups within a theoretical framework 
(Anisotropic Chromo Dynamics, ACD) that attributes to the fundamental quanta of QCD (Quarks 
and Gluons) a completely hadronic behaviour except, of course, for the Quark "point-like" coupling 
to the electroweak fields. 

The remarkable success of our calculation, free of adjustable parameters, highlights the fallacy of 
considering "short-distance" physics as autonomous (indeed orthogonal) from the long-distance, 
colour confining hadronic physics. 
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1 Introduction 



One of the most relevant and surprising aspects of the physics produced by HERA in the last four 
years is the discovery of a large number of events whose final states are characterised by large 
rapidity gaps , of the kind one observes in diffractive events in purely hadronic interactions . 
This occurrence is all the more surprising since "Diffraction" , or in general the physics of the 
Pomeron, is considered to be the genuine manifestation of the dynamics of QCD at large distances, 
or at small transverse momenta, where perturbation theory (PQCD) has absolutely nothing to say, 
the main actors being the Hadrons, the permanent prisons of Quarks and Gluons. On the other 
hand Deep Inelastic Scattering (DIS) is generally considered the true realm of QCD at short 
space-time distances where Asymptotic Freedom (AF) should allow us, should it not?, to compute 
the basic inclusive cross-section through the simple Feynman rules of PQCD. 
In hindsight a number of PQCD mechanisms Q have been invoked in order to save the general 
expectations based on AF, but the puzzle that stands in the way of the natural philosopher is to 
understand why a typical hadronic behaviour, completely extraneous to PQCD, can indeed be 
mimicked in a world thoroughly governed by this simple (unconfined) realization of QCD. 
Almost a quarter of a century ago one of us (CP.) found that the fundamental aspects of the just 
discovered Bjorken scaling behaviour could be easily and naturally accounted for by assuming that 
Quarks (and later Gluons) would behave at any space-time scale like hadronic particles [H, with 
their Regge behaviours that were known since long to well describe high energy interactions 
Naturally in the Massive Quark Model (MQM) ^ no attempt was made to derive "Reggeism" from 
any basic Lagrangian, the only firm Ansatz being that Regge behaviour, and in particular the 
Pomeron and Diffraction, should be properties of the Quarks as well. It should be perfectly clear 
that such an attitude was then (and is now) completely at variance with the beliefs and the 
expectations of the vast majority of the particle physicists; however it is a fact that the recent 
observations at HERA do agree with those distant attempts, while they have posed and continue to 
pose some non-trivial problems to the theorists of PQCD. 

While along the paths indicated by AF PQCD has developed into an extremely powerful and flexible 
means to deal with DIS physics, the MQM has gone through the much more uncertain steps that 
have led it to Anisotropic Chromo Dynamics (ACD) §, the "effective" Lagrangian which governs 
the dynamics of QCD over a peculiar confining vacuum (the Chromo Magnetic Liquid, CML |^) 
that has definite chances to well approximate the true QCD vacuum. Recently one has been able to 
find within ACD a simple and satisfactory approach to the Pomeron and Diffraction that 
embodies in a physically transparent way the pioneering approach of Low and Nussinov ||l0|| . 



2 Diffractive Deep Inelastic Scattering 

The basic idea is that when the Quarks of two Hadrons collide at high energy (see Fig.l) and interact 
by the fundamental colour-exchange potential, the Hadrons (the colour-singlets) of the final state 
comprise two different (in our case) QQ-pairs whose masses are of the order of the CM energies, 
which then fill the final states with the products of their decay . A detailed calculation shows 
that the diagrams in Fig.l correctly reproduce the main characteristic of the Pomeron. As for single 
Diffraction in DIS the relevant diagrams are reported in Fig. 2, where it should be noticed that the 
deep inelastic photon couples to the Quarks through a "point-like wave function" . 
The virtual photon-proton diffractive cross section is given by a tedious but straightforward 
calculation of the Feynman-like diagrams of Fig. 2, whose only "unorthodox" elements are: 

(i) the proton wave-function that projects the initial hadronic state (the proton) onto its 

constituent quark state which in the high energy limit (W^.p ^ oo,XBj ^0) can be taken as: 

// (2n)^2E 
d'md'p2^ (2^)92i^,2i;l2i;3 ^^^"^'^"^^ \P1,P2,P-P1-P2 > (1) 
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with the normalisation 

|dVlrfV2 |0(P1,P2)|'-1, (2) 

and the factorisation property {pi — Xip +p±i), \p±i\ <C \p\ 

(/) = 4)l{xi,X2) ■ (t>T{P±l,P±2), (3) 

which aUows the complete separation of the trivial longitudinal dynamics from the transverse 
one. We parameterise 0t(pi±,_P2±) as a Gaussian distribution 

(I)t{pi±,P2±) = Vl2- exp[-b{pl^ +pi^ + +pj_2)^)] (4) 

TT 

with b = 17 GeV~^, as determined by the slope of the forward peak in the pp elastic cross 
section. 

(ii) the diffractive final state comprises two Fire String (FS) states whose decay products are 
the Hadrons of the final states. Please note that in the inclusive cross-section Quark- Hadron 
"duality" O] allows us to integrate over the quark degrees of freedom only; 



(iii) the ACD gluon propagator 0| , whose form can be taken as 

, 3a^-\q^\ 
ia^ + \q^\r 



giving rise to an interquark effective potential V{r) = y?re^°'^ , where the two parameters /x, 
the effective string tension, and the "longitudinal size" of the "chromomagnetic needle" 
are fixed by the analysis of the hadronic spectrum iQ as {nig denotes the gluon mass) 

II = 0.475 GeV, a = 2V2mg ^ 0.40 GeV. (6) 
The result of the calculation can be expressed in terms of a diffractive structure function F^: 



drjpdt 



{xBj,Q ,ViP,t) (7) 



- XBj) d-qpdt 



where {xp is the Feynman variable of the final state proton, Mx is the invariant hadronic mass, t is 
the invariant momentum transferred to the final state proton). is given by the expression 

2 __J(i)2j(^^^.^Q2)(^_^^^.)^ (9) 



drjpdt rjjp 

where the normalisation contains the couplings, the colour factor, the sum over active quark flavour 
and the other numerical factors, 
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N^yQ)CF-(^) 1.77- 10-2 ggy-2 (^0) 

f 

I{t)^ j d^wV{v3-XQ)V{w + xo){e-^'-'^° -e-^""^) (11) 

with 3^0 = 2^ transverse momentum of the final state proton, and 5 = ba^ ~ 3. Furthermore 

r W\ 

J{x,Q')^ d^^d^zz^\n^V{z)^[cj,^{^,Q)-cj>^{^J)] (12) 
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where the "photon wave-function" (j)-y{^,z) is given by {A — ^): 



cl^,{tz)= I dr + ^ (13) 

' [r{l-r)A + eMl-r)A + {^ + zr] 







A few observations on Eqs. (5) to (9) are in order: 

(a) the energy dependence is logarithmic, thus no "unitarization" is necessary to obey the 
Froissart bound p5[ ; 

(b) in the hniited kinematic range now accessible to the experiments such logarithmic behaviour 
can well mimic a power law, which is favoured by the experimental fits; 

(c) the factorisation of the diffractive structure function as 

F,''{vjp,P,Q^) ^ ^rjjp) ■ Ffi(3,Q') (14) 

is broken due to the i-integration of equation (11) and to the contribution of diagrams 
subleading respect to that of Fig. 2, however, as one can see from Fig. ||, in the experimentally 
accessible region our result is effectively factorized; with factorizable form the F^{(i, Q^) can 
be associated with Pomeron structure function, as it is discussed in the pioneer work [ p^ . 

(d) when /3 ^ 0, corresponding to very large diffractive masses, the typical HERA energies 
(W^.p ~ 100 GeV) imply that rjp becomes rather sizeable, thus requiring the introduction of 
sub-leading Regge trajectories, like w,/2, ■■•,7r, as suggested by purely hadronic data \ YI_, D, 



that we have neglected in our calculation. Since such contributions can alter the form of the 
flux function $(r/jp), the extraction of {(3, Q^) can be more safely performed by integrating 
F^ in each bin of the (/3, Q^)-plane, rather than fitting the F^-function with a universal 
power law over a fixed rjp interval which is then used in the determination of Ff, as is 
commonly done by the HERA Collaborations |l9| . 



3 Comparison with experimental data 

Our calculation of the triple differential structure function F^ {rjp, l3,Q^) is compared with the 
experimental data from the ZEUS Collaboration in Fig. || and from the HI Collaboration |l^ in 
Fig. 1^. One notes good agreement with observations for all accessible values of P,Q^ and rjp, leading 
to overall statistical value x^/d.f. = 67.4/96 being compared with data. In order to see the 
consistency of both HI and ZEUS data we performed a new power-law fit to the combined set of 
data (shown in double-log scale in Fig. ^ and Fig. ^by straight line dependence) , which yielded the 
universal value of the Pomeron intercept ap{0) = 1.069 ± 0.030 {x'^/d.L = 33.33/47 [HI], 4.05/23 
[ZEUS] assuming full errors in the fit as statistical and systematic errors added in quadrature). 
Within experimental errors, this value of ap{0) is consistent with a soft supercritical Pomeron |^o|] 
and is correctly reproduced by our model. 

Our prediction for double differential structure functions F^{f], Q^) and F^{xBj, Q^) is reported for 
different values of in Fig. ^ and compared with data points from both ZEUS and HI [^8[ |l9|, pT| . 
Again good agreement is to be noted. In particular a smooth Q^-dependence is naturally reproduced 
by the model. To calculate the /3-distributions, the underlying ?7jp-integration is performed for two 
different intervals defined by HERA Collaborations [|l8|, |l^, and we find that the discrepancies in 
the values of the Pomeron structure function measured by HI and ZEUS does depend on such 
difference. Indeed, as one can see in Fig. ^ our model favours a small /3-dependence in the Pomeron 
structure function. The present (1993) data are statistically compatible with both the model of 
Ref. [|l8j and our predicted fiat /3-dependence. With more precise data it will certainly be possible to 
distinguish between the two approaches. But in this case we must be more careful in data analysis 
because the unfolding procedure (to calculate the acceptance and compensate kinematic cuts effects) 
based on Monte Carlo simulations [|l^, |l^, pl| can give undesirable correlations between the model 
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to be tested and the corrected data |18, since /3-distribution affects, in particular, the mass 
spectrum of hadronic final states. Clearly, a better comparison of our model could be made by 
processing with it the raw data. 

Finally, our expectation for the total contribution of Single Diffractive Dissociation (SDD) process is 
given in the last four pictures in Fig. in terms of scaling variable x bj and for different values of 
Q^. The relevant upper limit for the ?7jp-integration must be chosen about 0.15, in accordance with 
pure hadronic data ^. Note, that the decrease of to zero at XBj — 10^^ is due to the special 
kinematic and experimental conditions at HERA. When these limitations are removed it is seen that 
the real contribution of SDD to DIS is closer to 20% than to 10%, as claimed by the HERA groups. 
This discrepancy is due to their procedure for selecting large rapidity gap events and to the far from 
ideal acceptance of their setup for diffractive physics. 



4 Conclusion 

To conclude this letter, we may state with confidence that a sound, realistic and conceptually simple 
explanation has emerged from ACD of the complex features of Diffraction in DIS discovered by 
HERA. The main message of this paper, we believe, is the further proof of the correctness of the 
idea that deep inelastic phenomena fully exhibit the basic features of high energy hadronic 
dynamics, albeit at the level of the fundamental quarks instead of the composite hadrons. On the 
other hand the hadronic nature of the quarks that is corroborated by this work is, to our mind, a 
further stumbling block on the road of PQCD, which looks at the deep inelastic world as completely 
different and separated (by the scale Aqc-d) from the world of long-distance (low pt) phenomena. 
Finally we notice that the same ideas can be employed to compute the non-diffractive low XBj part 
of the structure function, which shall be reported in a future publication. 
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Captions 



Fig. 1 

The amplitude dominating the total cross-section for meson-meson scattering at high energy (all the 
diagrams obtained by permutations of the quark-gluon vertices should be summed). 

Fig. 2 

The amplitude dominating the single diffractive cross-section for DIS at high energy W^^-p and large 
diffracted mass (the so called Triple- Regge regime [12]). 

Fig. 3 

The triple difFcrcntial structure function of ZEUS is compared with our calculation (broken 
line), first passing through all the experimental cuts (dotted line) and then relaxing the cuts, but in 
the integration region defined by ZEUS (full line) and with a power-law fit (dashed line). The 
almost negligible contribution of non-factorisable diagrams is also shown (dot-dashed line). Note 
that the 15% contribution due to double diffraction has been subtracted from data [18]. 

Fig. 4 

The triple differential structure function of HI is compared with our calculation (broken line), 
first passing through all the experimental cuts (dotted line) and then relaxing the cuts, but in the 
integration region defined by HI (full line) and with a power-law fit (dashed line). Note that the 
15% contribution due to double diffraction has been subtracted from data [18]. 

Fig. 5 

The double differential structure function F^ of ZEUS and HI in terms of (3 and Xbj, predicted by 
the model (full line) and compared with the experimental data. The black points on /3-distribution 
correspond to ap{0) = 1.069, obtained in our power-law fit to combined set of HERA data. The 
black and white points on -distribution are observations using two different analysis methods 
[21], the vertical lines are pure kinematic cuts. Note that the 15% contribution due to double 
diffraction has been subtracted from data [18]. The last four pictures show the expected contribution 
of single diffraction dissociation, integrated with the relaxed upper limit on rip (see text above). 
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Figure 1 : The amplitude dominating the total cross-section for meson-meson scattering at high energy 
(all the diagrams obtained by permutations of the quark-gluon vertices should be summed) . 




► FS 

QHD - Cut 



Figure 2: The amplitude dominating the single diffractive cross-section for DIS at high energy W^'p 
and large diffracted mass (the so called Triple- Regge regime [12]). 



The Proton Diffractive Structure Function F2°(?7p,/3,Q^) I ZEUS Dote 




Figure 3: The triple differential structure function F2 of ZEUS is compared with our calculation 
(broken line), first passing through all the experimental cuts (dotted line) and then relaxing the cuts, 
but in the integration region defined by ZEUS (full line) and with a power-law fit (dashed line). The 
almost negligible contribution of non-factorisable diagrams is also shown (dot-dashed line). Note that 
the 15% contribution due to double diffraction has been subtracted from data [18]. 
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The Proton Diffroctive Structure Function F2°(?7p,/S,Q^) 1 HI Data 




Figure 4: The triple differential structure function of HI is compared with our calculation (broken 
line), first passing through all the experimental cuts (dotted line) and then relaxing the cuts, but in 
the integration region defined by HI (full line) and with a power-law fit (dashed line). Note that the 
15% contribution due to double diffraction has been subtracted from data [18]. 
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The Proton Diffractive Structure Functions F2°(/3,Q^) , F2°(x,Q^) 




Figure 5: The double differential structure function F2 of ZEUS and HI in terms of /? and XBj, 
predicted by the model (full line) and compared with the experimental data. The black points on 
/3-distribution correspond to aip(O) = 1.069, obtained in our power-law fit to combined set of HERA 
data. The black and white points on -distribution are observations using two different analysis 
methods [21], the vertical lines are pure kinematic cuts. Note that the 15% contribution due to double 
diffraction has been subtracted from data [18]. The last four pictures show the expected contribution 
of single diffraction dissociation, integrated with the relaxed upper limit on tjjp (see text above). 
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